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Abstract 

Since b quark mass is not asymptotically large, chirally enhanced 
corrections which arise from twist-3 wave functions may be important 
in B decays. We thus evaluate the hadronic matrix elements with the 
final light pseudoscalar mesons described by leading twist and twist-3 
distribution amplitudes. We find that chirally enhanced corrections can 
be included consistently in the framework of QCD factorization only if 
the twist-3 distribution amplitudes are symmetric. We then give explicit 
expressions of for i? — > vrvr at the next-to-leading order of including 
chirally enhanced corrections. We also briefly discuss the divergence 
appeared in the hard spectator contributions. 
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Non-leptonic two-body B decays are crucial for us to extract CKM matrix elements 
and uncover the origin of CP violations. Experimentally, with the running of B factories, 
there will accumulate a great amount of data on various B decay channels. Theoretically, 
however, how to extract CKM matrix elements from non-leptonic B rare decays with model- 
independence is still an open question due to the complexity of strong interaction. In the 
following, we first give a theoretical sketch on non-leptonic B decays. 

It is well known that the amplitude for the decay B P1-P2 can be expressed as 

A{B ^ P1P2) ^Y.^^C,{fi){P,P2m^l)\B), (1) 

i 

where Aj is a CKM factor, is a Wilson coefficient which incorporates short distance 

contributions from strong interactions and therefore is computable by making use of oper- 
ator product expansion and renormalization group equations, {PiP2\Qi{fi)\B) is a hadronic 
matrix element. Obviously, if we want to extract CKM factor from B decays, the hadronic 
matrix elements should be evaluated reliably. However, due to our ignorance on hadroniza- 
tion, it would be a great challenge to calculate these hadronic matrix elements reliably from 
first principles. A commonly used approximation is naive factorization assumption, which 
is based on Bjorken's color transparency argument |^: b quark decays and transfers a large 
momentum to final light quarks, in which two fast-moving, nearly collinear final quarks with 
appropriate color can be viewed as a small color dipole which will not significantly interact 
with the soft gluons and finally form an emitted meson. Then we have: 

(PlP2|g.(/i)|5) = (Pi|Ji|0)(P2|J2|5), (2) 

where Pi labels the emitted meson and P2 labels another light meson which absorbs the 
spectator quark from B meson. This approximation completely ignores non-factorizable 
contributions which connect the emitted meson to the spectator system and expresses the 
hadronic matrix elements in terms of meson decay constants and form factors. Since de- 
cay constants and form factors can be, at least in principle, well determined from other 
experiments, the branching ratios of non-leptonic B decays are obtained under this assump- 
tion. The main deficiency of this approximation is that non-factorizable contributions are 
completely missing. In consequence, the hadronic matrix elements lose their scheme- and 
scale-dependence. Noting that Wilson coefficients are scheme- and scale-dependent, the 
corresponding decay width will also depend on renormalization scheme and scale which is 
unphysical. This is a clear indication that non-factorizable contributions, which amount 
to final-state rescattering and strong interaction phase shift, are important. Several gen- 
eralizations of naive factorization assumption have been proposed to phenomenologically 
parameterize non-factorizable contributions. Since this kind of parameterization has no re- 
lation to, and therefore does not gain any information from, the underlying QCD dynamics, 
the resulting predictions on B decays are still mo del- dependent. 

In ref Beneke, Buchalla, Neubert and Sachrajda proposed a promising QCD factor- 
ization method: The hadronic matrix elements {PiP2\Qi{n)\B) contain two distinct scale: 
one is a large scale fi = 0{mb), the other is Aqcd which is the scale of hadronization. In 
the heavy quark limit, they show that the short distance contributions which are related 
to the large scale = 0{mb) can be, at least at one-loop order, separated from the long 



distance effects and are tlius calculable. Furthermore, the long distance effects can be pa- 
rameterized by light-cone distribution amplitudes and non-perturbative form factors. Thus, 
the factorization formula can be explicitly expressed as: [§,0 



1 X 



where and $Pi(a;) are the leading-twist light-cone distribution amplitudes of B and 

the final light mesons respectively, T/'^^ denote hard-scattering kernels which are calculable 
order by order in perturbative theory. This formula holds for the case that the emitted 
meson Pi is a light meson or an onia of two heavy quarks P,P|,pi| no matter whether 

P2 is a heavy or light meson. But in this article we will focus on the case that B decays to 
two light pseudoscalar mesons. 

In ref [|],|[, the authors pointed out that the equality sign of eq. (3) is valid only in the 
heavy quark limit. So if the heavy quark limit is an adequate approximation for B meson, 
or in another word, if power corrections in 1/mf, can be safely neglected, then everything is 
perfect. At the zero order of a^, it can reproduce "naive factorization", at the higher order of 
Os, the corrections can be systematically calculated in Perturbative QCD which will restore 
the scheme- and scale- dependence for the hadronic matrix elements. Therefore, the decay 
amplitudes of B meson can be reliably evaluated from first principles, and the necessary 
inputs are heavy-to-light form factors and light-cone distribution amplitudes. But in the 
real world, bottom quark mass is not asymptotically large(but about 4.8 GeV), therefore it 
may be necessary to consider power corrections in 1/mfe. Unfortunately there are a variety 
of sources which may contribute to power corrections in l/rrib, examples are higher twist 
distribution amplitudes, hard spectator interaction and transverse momenta of quarks in 
the light meson. Furthermore, there is no known systematic way to evaluate these power 
corrections for exclusive decays. Though naively, it is expected that power corrections 
may be neglected because AQcn/^b — 1/15 is a small number, power suppression may 
numerically fail in some cases. An obvious and possibly the most important case is chirally 
enhanced power corrections. As pointed out in ref 0], numerically the enhanced factor 

~ mhinT+ma) — ^-^^ which makes the power suppression completely fail. This parameter 
is multiplied by Og and ag, where is very important numerically in penguin-dominated 
B decays. So an evaluation of the hadronic matrix elements including chirally enhanced 
corrections may be phenomeno logically or numerically important. In the following, we will 
examine this problem in some details. 

Chirally enhanced corrections arise from twist-3 light-cone distribution amplitudes, gen- 
erally called ^p{x) and ^a{x). For light pseudoscalar mesons, they are defined as ^ 

{P{p my)^l5q{xm = fpfiP t du e^("P'-^+^^'-^)$,(«), (4) 

JO 

{P{p'myy,ul,q{xm = tfpMp',^u-pW)J^' du (5) 

where /ip = ^^j^^^ , z = y — x, mi and m2 are the corresponding current quark masses. If 
we want to generalize QCD factorization method to include chirally enhanced corrections 



consistently, we should describe the emitted light meson with leading twist-2 and twist-3 
distribution amplitudes P]: 
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A technical proof of factorization requires that the hard scattering kernels in Eq.(3) are 
infrared finite. Authors of Ref ||^,^ have shown it explicitly with leading twist distribution 
amplitudes. Then a basic and perhaps a difficult task for us is to show the infrared finiteness 
of the hard-scattering kernels using twist-3 distribution amplitudes after summing over the 
four vertex correction diagrams (Fig. l(a)-(d)). 

The start point for B decays is |Ai?| = 1 effective Hamiltonian 



Gp 
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+ h.c, (7) 



where Vg = VgbV*^(foT b ^ d transition) or Vg = V^fel/* (for b ^ s transition) and Cj(/i) are 
Wilson coefficients which have been evaluated to next-to-leading order approximation. The 
four-quark operators Qi are 

= {Uaba)v-A{qt3Ui3)v-A Q\ = {Caba)v-A{qf3Ci3)v-A 

Q2 = {Uabf3)v-A{qf3Ua)v-A Q\ = {Cabf3)v-A{q(}Ca)v-A 

Q3 = {qabaW-AEiq'fsq'fsh-A Qi = iqpba)v-AEiqaq'(})v-A 

q' q' 

Q5 = {qaba)v~Aj2{qi3q'f3)v+A Qg = {q/3ba)v-Aj2{qaq'l3)v+A (8) 

q' q' 

Qi = U^aba)v-AEegiq,3q'p)v+A Qs = U^(3ba)v-AEeg'{q^q'/^)v+A 

Q9 = liqaba)v-Aj2eq'iq/3qi3)v-A QlO = I{qi3ba)v-Aj2eg'{qaqi3)v-A 
q' q' 

and 

<577 = ^rnbqaCT^"'{l + -i5)bo.F^y, Qsg = ■^mbqa(y^"'tl.bpG'' {q = d or s). (9) 

With these effective operators, B P1P2 decay amplitudes in QCD factorization can be 
written as: 

A{B ^ P1P2) = % E E vX{PiP2mB)F, (10) 

V p=w,c j=l,10 

where Vp is CKM factor, {PiP2\Qi\B)p is the factorized matrix element. We will calculate 
QCD coefficients af and show explicitly that they are infrared finite. 

Infrared divergences exist in vertex correction diagrams (Fig.l(a)-(d)), so let us first 
consider these diagrams. For {V — A) 0(\^ — A) and [S + P) 0(5 — P) operators, twist-3 



distribution amplitudes make no contribution because of their Lorentz structures. Therefore, 
QCD coefficients af (except for Og and Og) are nearly as same as those obtained in Ref 
where leading twist distribution amplitudes are considered. The only difference is 
hard-spectator term (Fig.l(g)-(h)) which have been shown in Ref pJ]JT^ , we will discuss it 
later. As to (l^ + A) — A) operator, there are some subtleties in regularizing the infrared 
divergences. If we use dimension regularization, the infrared finiteness will not hold after 
summing over those four vertex correction diagrams. That is because wave functions are 
defined in 4-dimensions, it may be unconsistent to naively extend its usage to d-dimensions. 
Thus we assign a virtual mass to the gluon propagator and regularize the infrared integrals 
in four dimensions. For the twist-3 distribution amphtudes $p(x), the calculations are 
performed in momentum space. Then it is straightforward to verify that the vertex correction 
contributions of {V + A) 0(1^ — A) operator to {S + P) 0(5' — P) are infrared finite: 



V - - — — 
Att N 



/ dx^p{x){iTi\og- + \og- - Li2{ — ) + Li2(--) + 6}, (11) 



where x = 1 — x and Li2{x) is dilogarithm function. On the other hand, when considering 
we have to do the calculations in coordinate space according to Eq.(5). For example, 
let us consider Fig. 2. In coordinate space, we have: 



A = J d'x,d'x2d'x, w,e^«'-^3 ^,(1 + ^5) j ^^e-'^-^^^-..) [-^g^^-{T^ 
'^''^''-a,.l,q'{xs - x^r I dve'^^'^-^^''^-^^^4^5p. - 75) 
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(27r)4 kl-ml ' ^^'"^ ^'^"^ " J (27r)4A;2-mr 



^-^^qf'ixs - X2y^^6f3^e''^'"'-'"'+'"'''^^ g-ifci.(x3-xi) e-'P-'^^ba- H^^{k,kuk2) 

47V r- 6 



fca, (12) 



k2=k—vq 



where H^i,(k, ki, k2) contains the Lorentz structure and propagators of the hard scattering 
kernels: 



H,Ak, k,, k2) = 7a(1 + 75) S[-^^7s7"(T")p.] 7'(1 - 75) 

/Co 



"^2 



i(^i + mfe) . -i . . 

[-^^^7«(T ),„] ^^V- (13) 

After a lengthy derivation, we can regularize the infrared divergences with a gluon virtual 
mass mg-. 

Fig.l(a) ~ — ^ ^ { ^ ^ + 21og(—f) log /i — 4 log flog /i + log /i + finite terms} , (14) 
Atx Nc V 2 



Fig. 1(b) ~ '^^^^ r A'- _|_ 2 iog(— ti) log /i — 4 logf log /i + log /i + finite terms} , (15) 

4:71 Nc V 2 

Fig. 1(c) ~ ^-^—^^^{log^ /i — 2 log(—t') log /i + 3 log /X + finite terms} . (16) 
Fig. 1(d) ~ —^^^—^i^^{log^/i — 2 log(—'u) log + Slog + finite terms} , (17) 

where n = m'^/ml. From the above equations, it is observed that, in the case of $a- distri- 
bution amplitudes, the terms with infrared divergence in vertex correction diagrams can not 
cancel unless ^a{v) is a symmetric function: $o-(f) = ^a{v). This is an unexpected result, 
which means QCD factorization is violated for asymmetric twist-3 light-cone distribution 
amplitudes. This indicates that chirally enhanced corrections can be included consistently 
in the framework of QCD factorization only when twist-3 light-cone distribution amplitudes 
are symmetric. Therefore, in the following, we will implicitly assume a symmetric twist-3 
light-cone distribution amplitude for light pseudoscalar mesons. It is then straightforward 
to show that vertex corrections of {V + A) ^ {V — A) operator are completely canceled after 
summing over four diagrams in the case of $cr distribution amplitude. 

For penguin contractions (Fig.l(e)-(f)) and hard spectator diagrams (Fig.l(g)-(h)), we 
shall also do the calculations in coordinate space when $o-('y) is included. When treating 
penguin contractions, it should be careful that Fig. 1(e) contains two kinds of topology, which 
is displayed in Fig. 3. They are equivalent in 4 dimensions according to Fierz relations. 
However, since penguin corrections contain ultraviolet divergences, we must do calculations 
in d dimensions where these two kinds of topology are not equivalent [1^. We did not 



notice it and therefore obtained a wrong term — ^C*4 in the expression of 04 in 0. We also 
obtained a wrong term (C3-|-C4/A^)/3 and missed a term of {CA + C3/N) in the expression 
of Oj'q in for the same reason. 

Then as an illustration, the explicit expressions of (i = 1 to 10) for B ^ titt (using 
symmetric light-cone distribution amplitudes of the pion) are obtained. But it is easy to 
generalize these formulas to the case that the final states are other light pseudoscalars. We 
now list af for S — vrvr as follows: 



-Tlf{^4'°^- + - |) + (C3 - log ^ + G(0) + G(l) - 1) 

+ E iC, + C, + le^Cs + le^C,o)dlog-^ + Gis,)) + GsCsG]: (21) 

q=u,a,s,c,b " 

= ^= + § + - 12). (22) 
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+ E (^3 - ^)((1 + log ^ - - + G\s,) + G'^l.,)) 



+ E (<^4 + + ^eqCg + ^eqCio) ((1 + \a,) log - ^ - + G\s,) + G'^ls, 

q=u,d,s,c,b 

+{l + A,)Csa} , (23) 
-^-^r + ^+'^^C,(-F-ll), (24) 

{(C. + f )((! + f ^„) log|^ - ^ - + + G'(..)) 

+(C. + f ) E ja,((l + Ia,) log £^ - ^ - + G'(.,) + G'(s,)) 

q=a,b " 

+(C3 + ^ + a + §) E Je,((l + HA,)logA-±-l^„ + G'(.,)+G'(.,)) 

q=u,a,s,c,b ^ o 

+(| + ^^.)C7.} , (25) 

a9 = C^9 + ^ + ^^Cioi^, (26) 

a?o = ^10 + ^ + ^^^^9F -^[iC. + ^)(- log - + G{s,) - -) 

+(^4 + §) E Uihog-^ + G{s,)-l) 
^ g=d,b^ 3 mft 3 

+iC3 + ^ + C, + ^) ^e,(^log^ + G(.,)) + ^G8C,J. (27) 

q=u,d,s,c,b " 

Here = 3 is the number of color, Cp = ^2n^ is the factor of color, Sq = rrt^/ml and we 
define the other symbols in the above expressions as: 

F = -121n-^ - 18 + /^ + /^^ (28) 
rUb 

1 1 
f = J dxg{x)^x), Gs = J dxGs{x)^x), (29) 



1 

G(s)^J dx G{s,x)^x), (30) 



1 

G'{s) = J dx G'{s,x)%{x), (31) 



dx G (s,x)—— r, — [ dx — — r 

6 1 — a: J 61 — X 

^ ^ 



(32) 



here $(a:)($p(x),$o-(2:)) is leading twist (twist-3) wave function of the emitted pion, and the 
hard-scattering functions are 



1 — 2x 

g(x) = 3 Inx — SiiT, Gs(x) 

1 — X 



1 — x' 



1 

G{s, x) = —4 J du u{l — u) ln(s — u{l — u){l — x) — it) 



(33) 
(34) 



G"(s,x) 



-3 j du u{l — u) ln(s — u{l — u){l — x) — ie), 



(35) 



1 

G''{s, x) = -2 J du u{l - u) ln(s - u{l - u)(l 



X] 



+ 



J. 

/ 



du 



u^{l-uY{\-x) 
- u{\ — %L){\ — x) — ie 



(36) 



The contributions from the hard spectator scattering (Fig.l(g)-(h)) are reduced into the 
factor /^^. We take the wave function of B meson as 75 (-fe ~ Then an exphcit 

calculations show that twist-3 distribution amplitudes of the emitted pion make no contri- 
butions to /^^. It means that there is no hard spectator contributions for ag and ag. For 
other QCD coefficients af , we have: 



/ 



47r^ 



B 



N Ff^^(0)m' 



B 



d^ 



dx 



X 



dy 



$(?/) , 2/i^ $,(?/) 



l-y Ms 6(1 -y)^ 



(37) 



Here $(a;) is leading twist distribution amplitude of the emitted pion, $(?/)($o-(?/)) is twist- 
2 (twist-3) distribution amplitudes of the recoiled pion. This formula is consistent with the 
result of Ref. 0. 

In the above expressions of , Og and can now be evaluated to next-to-leading order of 
as, which significantly reduce their scale-dependence. As to other QCD coefficients af , there 
contains a divergent integral in hard spectator term /^^. In the next paragraph, we will argue 
that this disturbing divergence may need further consideration. Here we simply assume that 
/ ^ ~ In (similar to what have been done in Ref |lll],12|, though our assumption here 

is certainly an oversimplification). We thus illustrate numerically the scale-dependence of 
af(7r7r) in Table. 1. Here we use the asymptotic distribution amplitudes 



$(x) = $^(x) = 6x(l 



X] 



and ^p{x) 



(38) 



and the input parameters are taken as follows: F-^'^(O) = 0.33, = 0.2 GeV , = 
133 MeV , the pole masses = 4.8 GeV , m^. = 1.4 GeV, the MS masses Witirnt) = 
170 GeV, mb{mb) = 4.4 GeV, m„(2 GeV) = 4.2 MeV, md{2 GeV) = 7.6 MeV and 
A^^Ij, = 225 MeV. 



We notice that the above approach of evaluating hard spectator contribution is naive. 
For instance, the scale of hard spectator contribution should be different from the vertex 
correction contribution. While it seems reasonable to take the scale n ~ (9(mb) for the 
vertex correction diagrams to avoid large logarithm a^log^, a natural choice of the scale 
of hard spectator contribution may be around 0{1 GeV) because the average momentum 
squared of the exchanged gluon is about 1 GeV^. Another disturbing feature of hard spec- 



tator contribution is that, as pointed out in ref IITIIJT4 , when including the contribution of 
$0-, there would appear divergent integral Jq dy^ even if the symmetric distribution ampli- 
tude is applied. This divergent integral implies that the dominant contribution comes from 
the end-point region, or in another word, it is dominated by soft gluon exchange. How- 



ever the transverse momentum may not be omitted in the end-point region [jl4[, if so, the 
corresponding divergent integral would then changed to: 

! iyl ^ ! iyi-kr^^. (39) 

J y J y^mi + kj. 

As an illustration, we do not consider the /c^ dependence of wave functions (though it is 
certainly not a good approximation), then the above integral is proportional to: 

dydk^ f dxdy 

— : — n TT oc / . (40) 

yc,m( + kj. J X + y 

The above integration converges now, furthermore it is not dominated by end-point con- 
tribution. This illustrates that the treatment of hard spectator diagrams may need further 
discussions. 

There exists " annihilation" contributions which may belong to chirally enhanced correc- 



tions. In Ref. |]I2[, the authors have discussed this topic and find that a divergent integral 
(/ if )^ will appear. We suspect that this divergence may disappear, similar to the hard 
spectator term, if the effect of transverse momenta can be included. It is also possible 
that "annihilation" contributions are really dominated by soft interactions and thus violate 
factorization. Due to its complexity, we do not include "annihilation" contributions in the 
expressions of af. 

In summary, to generalize QCD factorization method to include chirally enhanced correc- 
tions consistently, the final light mesons should be described with leading twist and twist-3 
distribution amplitudes. We demonstrate that the infrared finiteness of the hard scattering 
kernels can be obtained only if the twist-3 distribution amplitudes are symmetric. We then 
give explicit expressions of af at next-to-leading order of a.^ including chirally enhanced 
corrections. We also discuss briefly the disturbing hard spectator contributions. 
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TABLES 





^ = 5. (J Gel/ 




fi = 2.0 GeV 




uoemcients 


JNLU 


LU 


ATT 

JNLU 


T 

LU 


a" 


1.024 + O.Olzz 


1.017 


1.034 + 0.024z 


1.037 


of 


0.144 - 0.076i 


0.188 


0.123 - O.lOOi 


0.109 


as 


0.003 + 0.002i 


0.002 


0.004 + 0.004i 


0.004 


al 


-0.027 - 0.014i 


-0.029 


-0.029 - 0.017i 


-0.040 


a% 


-0.033 - O.OOTi 


-0.029 


-0.036 - 0.007? 


-0.040 


as 


-0.003 - 0.003i 


-0.005 


-0.002 - 0.005? 


-0.010 




-0.036 - 0.012i 


-0.033 


-0.037 -O.Ollz 


-0.040 




-0.039 - O.OOSi 


-0.033 


-0.040 - 0.004i 


-0.040 


a7 X 10^ 


11.9 + 2.8i 


13.8 


0.0 + 5.4i 


7.6 


r-^a^ X 10^ 


36.8 - lO.Qz 


36.8 


45.0 - 5.2? 


39.8 


r^a§ X 10^ 


35.0 - 6.2i 


36.8 


44.2 + 3.1? 


39.8 


ag X 10^ 


-936.1 - 13.4i 


-928.4 


-953.9 - 24.5? 


-957.3 


0^*0 X 10^ 


-81.8 + 58.8i 


-141.4 


-58.3 + 86.H 


-74.0 


afo X 10^ 


-85.2 + 63.5i 


-141.4 


-60.3 + 88.8? 


-74.0 



TABLE I. The QCD coefficients (vrvr) at NLO and LO for the renormahzation scales at 
it = 5 GeV and u = 2.5 GeV, where = — — ^ 



FIGURES 




(a) (b) (c) (d) 




(e) (f) (g) (h) 



FIG. 1. Order of ag corrections to hard-scattering kernels. The upward quark hnes represent 
the ejected quark pairs from b quark weak decays. 




FIG. 3. Two kinds of topology for penguin contractions. 



